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ABSTRACT 


Presented in this report are both the profile and bulk aerodynamic derived micrometeoro- 
logical data taken in the marine atmospheric surface layer at San Nicolas Island, California, dur- 
ing a five day experiment in November 1978. Thirty-six figures are presented to represent 
much of the data results in graphic form. A comparison of the profile and bulk derived stability 
parameters strongly indicated an escarpment phenomenon similar to what had been reported for 
the Risé, Denmark, site. Although an escarpment induced, height dependent horizontal 
acceleration of the wind field profile was observed, there was no apparent affect on either the 


temperature or humidity profile. A suitable wind field escarpment correction for future data 


taking was left to the next experiment, in which a better quality bulk water temperature meas- 


urement essential to the bulk comparison would be available. 








MICROMETEOROLOGICAL DATA REPORT FOR 
THE NOVEMBER 1978 ELECTRO-OPTICS METEOROLOGY 
(EOMET) EXPERIMENT AT SAN NICOLAS ISLAND, CALIFORNIA 


INTRODUCTION 


The purpose of the Naval Research Laboratory’s (NRL) San Nicolas Island meteorological 
experiments is to provide the micrometeorological measurements necessary to characicrize the 
marine atmospheric surface layer in which simultaneous infrared atmospheric transmission 
measurements are being performed. The original impetus for this work was provided by the 
Optical Signatures Program of the Naval Weapons Center and has subsequently been 
transferred to the auspices of the Electro-Optics Meteorology Program (EOMET) of the Naval 


Ocean Systems Center. 


This report contains the results of one hundred and thirty-nine data runs constituting 
sixty-nine and a half hours of data taken over a five day period between 13 and 17, November, 
1978. Since the previous experiment (CEWCOM-78) of May, 1978, several changes and 
modifications have been incorporated into the present experiment. They include: 1. The 
introduction of a bulk aerodynamic technique for estimating the stability and flux parameters. 
2. The presentation of many of the parameters as a function of time in graphic form to facili- 
tate utilization. 3. The detailed graphic presentation of the profile and bulk stability 
discrepancy in order to determine and evaluate the degree of an island escarpment effect on the 
data. 4. The presentation of the parameter calculations in a simple algebraic form utilizing a 
short hand english notation as well as a listing of the computer program in the Hewlett-Packard 


HPL language. 


The experiment data contained in the present report is available on digital magnetic flexi- 


ble disk (Hewlett-Packard HPL language formatted) upon request to the NRL Atmospheric 


Physics Branch. The complete text of the report (including experiment data and figures) is 
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available on microfiche and information on how to obtain copies is available from the NRL 


Library Documents Section. 


In addition to the data contained in this report, there are two other NRL November 1978 
EOMET experiment memorandum reports. They cover the Radon 222 data taken by R.E. Lar- 
son for determining open ocean versus coastal conditions, and the aerosol distribution data 


taken by R.K. Jeck. 








LISTING OF DATA BY RUN NUMBER 


The Run Number signifies the data run start time: year, month, day, hour, and minutes, in 
Pacific Standard Time. 


13 Nov. 1978 14Nov. 1978 15 Nov. 1978 16 Nov. 1978 17 Nov. 1978 
7811132135 7811140000 7811150330 7811160000 7811170000 
7811132200 7811140030 7811150400 7811160030 7811170030 
7811132230 7811140100 7811150430 7811160100 7811170100 
7811132300 7811140130 7811150500 7811160130 7811170130 
7811132330 7811140200 7811150530 7811160200 7811170200 
7811140230 7811150600 7811160230 7811170230 
7811140300 7811150630 7811160300 7811170300 
7811140330 7811150700 7811160330 7811170330 
7811140400 7811150730 7811160400 7811170400 
7811140430 7811151301 7811160430 7811170430 
7811140500 7811151330 7811160500 7811170500 
7811140530 7811151400 7811160530 7811170530 
7811140600 7811151430 7811160600 7811170600 
7811140630 7811151532 7811160630 7811170630 
7811140700 7811151600 7811160700 7811170700 
7811140730 7811151630 7811160805 7811170730 
7811140800 7811151700 7811160830 7811170800 
7811140930 7811151730 7811160900 7811170830 
7811141000 7811151800 7811160930 7811170900 
7811141030 7811151900 7811161000 7811170930 
7811141400 7811151930 7811161030 7811171000 
7811141430 7811152000 7811161100 7811171030 
7811141500 7811152030 7811161130 
7811141530 7811152100 7811161200 
7811141600 7811152130 7811161230 
7811141630 7811152200 7811161300 
7811141700 7811152230 7811161330 
7811141730 7811152300 7811161440 
7811141800 7811152330 7811161500 
7811141830 7811161530 
7811141900 7811161600 
7811141930 7811161630 
7811142000 7811161700 
7811142030 7811161730 
7811142100 7811161800 
7811142200 7811161830 
7811142300 7811161900 
7811161930 
7811162000 
7811162030 
7811162100 
7811162130 
7811162200 
7811162230 
7811162300 


7811162330 








EXPERIMENT PROCEDURE AND NOTES 


The basic system description and rationale for this experiment are contained in the data 
report! of the previous May 1978 experiment. Since that time, the following modifications to 
the system and additions to the calculation procedures were incorporated into the present 
experiment. For a detailed presentation of the updated calculation procedures, see the CAL- 


CULATIONS section of this report. 


1. In order to facilitate the comparison of the micrometeorological and atmospheric 
infrared transmission data, each data run was identified by a mutually employed ten digit Run 
Number which represented the data run start date and time in terms of Pacific Standard Time 
(PST). For example: Run Number 7811131430 represents a data run started at 1978, 


November, 13th, 14:30 PST. 


2. In order to facilitate the synchronization of the simultaneous micrometeorological and 
infrared transmission data, the micrometeorological data averaging period was set for a duration 


of 30 minutes with each data run initiated at the clock hour or half hour. 


3. Following the practice of other experiments? conducted on a similar beach or island 
type escarpment, this experiment represented the measurement height above the escarpment, 
Figure 1, as being equivalent to the measurement height above the ocean. For simplicity, the 
indicated profile measurement heights (Z1 and Z2) were defined as the average height of the air 
temperature, wind speed, and dew point sensors at the respective levels. The turbulence meas- 
urement heights (ZTUR1 and ZTUR2) were located 0.60 meters above the respective profile 


measurement heights. 


4. In order to enable a comparison in this experiment of the profile and the bulk aero- 
dynamic techniques for computing the various micrometeorological parameters, the bulk aero- 


dynamic equations of C. A. Friehe’* (with minor modifications) were added to the post- 


experiment data reduction procedure. 
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Figure |. The measurement height above the island escarpment (Z, ) is equivalent to the 
measurement height above the ocean (Z,,). 





5. In order to compensate the dew point electronics (excluding the sensor) for time 
dependent variations and the non-interchangability of the device’s component subsystems, high 
precision resistance type secondary dew point standards (calibrated at the time of each sensor’s 
primary NBS calibration) were employed in the field. This procedure was typically performed 
twice a day to yield field calibration corrections (DPIFCAL and DP2FCAL). These corrections 
were added to the respective dew point raw data voltage value before employing the curve fit 
equations to convert the data into engineering units. It was found that when the dew point 
electronics functioned properly, the electronics typically drifted an equivalent of ~0.02°C (or 


~1/3 the magnitude of the measurement accuracy) over a one week period. 


6. While conducting the above described field calibration procedure, it was observed that 


« 


the top level dew point device (DP1) experienced drifts equivalent to ~0.20°C (or ~3 times 


the magnitude of the measurement accuracy) over only a 6 hour period. In a later checkout of 
6 


the dew point electronics, it was found that the erratic drift was due to a failure of the device’s 
sensor control circuit. As part of the post-experiment data reduction, an attempt was made to 
correct the data by assuming that the nonsystematic drift of the system was linear as a function 
of time between the field calibrations. However, as can be seen from the corrected data con- 
tained in this report, the specific humidity values of the upper level often exceeded those of the 


lower level, leaving the accuracy of the top level dew point readings still in question. 


7. In order to minimize the affect of the questionable accuracy of the upper level dew 
point measurement in computing the profile derived stability parameters (and thus the profile 
flux and profile scaling parameters), the profile Richardson’s number stability parameter was 
defined in terms of the potential temperature partial derivative and the ambient virtual potential 


temperature. 


8. The bulk aerodynamic calculations employed for this experiment were based upon the 
parameter values at the standard 10 meter height. The various parameter values for the 10 
meter height were estimated from the profile measurements by assuming an ideal log-linear 
profile relationship and taking into account the fact that each type of sensor system was slightly 


offset from the indicated profile measurement height. 


9. Unfortunately, the bulk water temperature measurement for this experiment still relied 
on the "bucket and thermometer" technique. In the post-experiment data reduction procedure, 
the raw data voltage value for the bulk water temperature was manually altered by assuming a 
linearly interpolated value for the data runs which were conducted between the actual bucket 


measurements. See figure 3 on page 12 and ERRATA NOTE | on page 60. 


10. The MRI model 1022 S anemometer cup failure experienced during the previous 
experiment was corrected by a special electron beam welding modification made by the 
manufacturer. It had been found that the solder which held the previous stainless steel cups to 


the anemometer turned to a powdery consistency in as little as three days exposure to the 


7 











marine environment. This caused the cups to rotate out of correct vertical alignment and/or 
fall off of the anemometer. The data contained in the previous experiment data report! had 
been screened to remove any data affected by this problem. The modified anemometers 
employed in this experiment were calibrated at the National Bureau of Standards in January and 
March 1979. To facilitate a greater measurement accuracy at low wind speeds, the calibrations 
and the representative polynomial curve fit equations were divided into two discrete segments: 
0 to 5 m/sec and 5 to 20 m/sec. This reduced the combined (calibration and curve fit) error to 


less than +2% of the reading for wind speeds above 2 m/sec. 


11. The R. M. Young model 43404 air temperature solar radiation shield back flow prob- 
lem, observed during the previous experiment, was corrected by the installation of the more 
powerful model 43404A aspirator onto the existing radiation shields. It had been observed that 
under high wind speed conditions, the mouth of the shield containing the air temperature sen- 
sor acted as a Bernoulli effect, or Prandtl tube, causing the air flow through the aspirator t@ 
reverse. Thus, air heated by the aspirator motor w.. drawn past the air temperature sensor and 
gave rise to spurious air temperature readings. The previous experiment report! did not contain 
data affected by this problem. The new aspirators employed in the present experiment also 
incorporated a "correct flow" indicator lamp powered by an internally mounted air flow detector 


switch. 


12. The YSI-Sostman barometric pressure sensor had been damaged in shipment to the 
previous experiment. In preparation for the present experiment, it and the backup sensor sub- 
sequently employed in the previous experiment were both returned in September 1978 to the 


manufacture for repair and recalibration. 


13. The polynomial curve fit equations (employed for converting the raw data voltages 
into engineering units) for wind speed, dew point, and barometric pressure, were updated and 


simplified by utilizing a library computer program on the NRL DEC-10 computer. 








14. The correction of both the profile and bulk derived specific heat of «ir as a function 








of specific humidity was added to the calculations. 


15. The correction of both the profile and bulk derived latent heat of water vaporization 








as a function of air temperature was added to the calculations. 






16. The computation of both the profile and bulk derived drag coefficient for the 10 


meter height was added to the calculations. 








17. The computation of the profile derived Z/L (or &) stability parameter for the 10 


meter height was added to the calculations to facilitate a comparison of the profile and bulk 








aerodynamic techniques at a standard height. 


18. Because the profile derived roughness length is a function of the profile Richardson 






Stability measurement, it is most valid at the geometric mean height. Thus, the profile rough- 






ness length computation was modified to be a function of the wind speed at the geometric 






mean height. 







19. The computation of the buoyancy flux was eliminated and the computation of the 








total heat flux budget (latent heat + sensible heat — sky and solar radiation) was substituted. 






20. The computation of vapor pressure in terms of torr was eliminated. 





21. Although no microthermal turbulence temperature structure (C,2) data were subse- 






quently taken during the present experiment, the capacity to compute the optical refractive 


index (C,?) from the turbulence data was added to the calculations. 







22. The computation of the tide table value was modified to be expressed in terms of 






meters and referenced to the mean sea level. 
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PARAMETER GRAPHICS AS A FUNCTION OF TIME 


Air temperature (ATB) 
Bulk water temperature (WTB) 


Potential air temperature minus bulk 
water temperature (PTWTB) 


Dew point (DPB) 

Relative humidity (RHB) 

Specific humidity (SHB) 

Absolute humidity (AHB) 

Vapor pressure (VPB) 

Wind speed (WSB) 

Bulk friction or scaling velocity (FRICVB) 
Bulk roughness length (RUFLB) 

Bulk drag coefficient (DRAGB) 
Barometric pressure 

Wind direction referenced to true North 
Sky and solar radiation 


Tide table referenced to mean sea level 


The following meteorological parameters were computed for the standard height of 10 


meters and plotted as a function of start time (Run Number) in PST. See pages 11-26. 


figure 2 


figure 3 


figure 4 
figure 5 
figure 6 
figure 7 
figure 8 
figure 9 
figure 10 
figure 11 
figure 12 
figure 13 
figure 14 
figure 15 


figure 16 


figure 17 
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PROFILE AND BULK STABILITY ANALYSIS: 


AN ESCARPMENT EFFECT 


The profile measurement technique is considered to be more desirable than the bulk aero- 
dynamic estimation technique for characterizing the atmosphere because it entails fewer 
assumptions and approximations. However, the much simpler bulk technique can play an 
important role in estimating the order of magnitude and sign of the stability and flux parame- 
ters. This experiment has brought to light a much larger than anticipated discrepancy between 
the profile and bulk derived stability values. What constitutes a reasonable agreement between 


the profile and bulk data is still, however, an unresolved issue. 


The accurate measurement of the profile Richardson stability parameter is the “linchpin" 
of the profile technique for computing the various flux parameters. Thus in this report, we 
confined our analysis to the apparent discrepancy in the Richardson stability values. To facili- 
tate this analysis, Figs. 18 through 36 (pp. 32-50) were prepared in an attempt to determine the 


most likely causes for the discrepancy. See ERRATA NOTE 3, page 60. 


Figures 18 and 19: the profile (RI) and bulk (RIB) Richardson stability as a function of 
time (Run Number). The profile derived stability value generally tracked the bulk value but 
was subjected to significantly larger magnitude perturbations from the neutral or near zero state. 
There were two major discrepancies in the stability sign which occurred on 15 November 
(3:30-7:30) and 17 November (9:00-10:30). Since the profile technique contains fewer sweep- 
ing approximations than the bulk, it was not unexpected that the profile values contained more 
reading to reading variation or fine structure about a six hour mean. However, this could not 
account for the typical order of magnitude greater accentuated features in both the positive and 


negative values of the profile stability data. 


Figure 20: the profile and bulk stability difference referenced to the bulk value 


((RIB — RI)/RIB) as a function of time (Run Number). If the value of +50% agreement were 
27 





arbitrarily chosen as the maximum difference of what could be considered a reasonable agree- 
\ ment between the profile and bulk values, then the two values were in agreement only 26% of 


the time. 


| Figures 21 and 22: the profile (RI) and bulk (RIB) stability as a function of the wind 
speed (WSB) at the 10 meter height. Above a wind speed of 7 m/sec, the two stability parame- 
ters displayed a similar systematic relationship to the wind speed. Below 7 m/sec, however, the 
profile stability data tended to divide into two diverging legs, one positive and one negative. 


Each leg was of greater relative magnitude than its single leg bulk counterpart. 


Figure 23: the correlation between the profile (RI) and bulk (RIB) stability values. There 


was little correlation between the two values, except under neutral or near zero conditions. 


Figure 24: the profile and bulk stability difference as a function of the profile potential 
temperature differential (PT1—PT2). The potential temperature differential (plotted with the 
sign reversed) determines the sign of the profile stability and to a lesser extent its magnitude. 
The typical maximum differential temperature error was +.006°C. No discernible pattern in 


the data was detected. 


Figure 25: the profile and bulk stability difference as a function of the profile dew point 


differential (DP1—DP2). The dew point differential (plotted with the sign reversed) played 





only a minor role in the profile stability, since for this experiment the humidity was included 
only in the ambient virtual potential temperature term and not in the profile temperature partial 
derivative. The typical maximum differential dew point error was +0.14°C. No discernible 


pattern in the data was detected. 


Figures 26 and 28: the profile and bulk stability difference as a function of the profile 
wind speed differential (WS1—WS2). The square of the wind speed differential appears in the 
denominator (1/(WS1 — WS2)2) of the profile stability Richardson number equation. Although 


it played no role in determining the stability sign, it played the dominant role in determining 
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the profile stability magnitude. In general, it was found that the smaller the differential, the 
larger the stability discrepancy. The largest stability discrepancies occurred when the wind 
speed differential was less than 0.5 m/sec. A typical maximum differential wind speed meas- 


urement error was +4% of the ambient wind speed. 


Figure 27: the profile and bulk stability difference as a function of the wind speed (WSB) 
at the 10 meter height. In general, it was found that the lower the ambient wind speed the 
larger the stability discrepancy. The largest stability discrepancies occurred when the wind 


speed was below 7 m/sec. 


Figure 29: the profile and bulk stability difference as a function of wind direction. 
Although the island escarpment upon which the readings were taken could be said to be very 
crudely symmetrical as a function of wind direction (to 90° either side of the Northwest), the 
path length of the wind field over the rocks between the escarpment and the water varied 


Significantly as a function of wind direction. No discernible pattern in the data was detected. 


Figure 30: the profile and bulk stability difference as a function of sky and solar radiation. 


It had been speculated that convective heating, particularly under low wind speed conditions, or’ 


reflected solar radiation from the island escarpment may have affected the temperature profile | 


measurements. No discernible pattern in the. data was detected. 


Figure 31: the profile and bulk stability difference as a function of tide height. The path 


length of the wind field over the rocks between the measurement escarpment and the water 


varied from zero to one hundred meters in length depending upon the tide. No discernible pat-“ 


tern in the data was detected. 


Figure 32: the profile and bulk stability difference as a function of the potential air tem- 
perature at 10 meters minus the bulk water temperature (PTWTB). No discernible pattern in 


the data was detected. 
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Figure 33, 34, and 35: the profile and bulk stability difference as a function of the bulk 
derived friction velocity (FRICVB), drag coefficient at 10 meters (DRAGB), and the roughness 
length (RUFLB). In general, it was found that when these three bulk parameters decreased, 
the stability discrepancy increased. However, this same relationship was observed with the 
ambient wind speed and the above patterns were due to the way in which these bulk parameters 


are "coupled" to the wind speed measurement. 


Figure 36: the profile and bulk stability differences as a function of the profile derived 
roughness length (RUFL). In general, the smaller the profile derived roughness length, the 
larger the stability discrepancy. The values ranged from one to twelve orders of magnitude (') 
smaller than 1 x 10-5 meters. Compare figures 35 and 36. Values ranging smaller than 1 
x 1075 meter are physically unrealistic and were a conspicuous telltale symptom of the root 
cause. The same type of results were observed at the Risé, Denmark, test site and reported by 


Panofsky and Peterson’. 


From the physically unrealistic telltale profile roughness length values reported by Panof- 
sky and Peterson®, Jensen and Peterson2 successfully accounted for the phenomenon by 
demonstrating an escarpment induced height dependent horizontal acceleration of the wind field 
profile. They further demonstrated that for the heights and distances such as those which exist 
at the San Nicolas Island test site, the lower level wind speed measurement would be 
accelerated more than the upper level measurement. This results in steeping the wind speed 
profile slope, decreasing the wind speed differential, and increasing the relative magnitude of 
the profile stability value independent of its sign. This explanation corresponded to the results 


demonstrated by our stability analysis. 


Although there is, as of yet, no suitable explanation for the two instances of profile stabil- 
ity sign reversal, there is no evidence in this experiment to indicate any detectable influence of 


the island on the temperature or humidity field profiles. 
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Due to the fact that an accurate and reliable bulk water temperature measurement is cru- 


cial to the bulk aerodynamic technique, and that such a system would not be available until the 


next field experiment (May 1979), it was decided to limit the scope of this report to a determi- 


nation and evaluation of the island escarpment problem. It would be left to the next experi- 


ment to provide a suitable correction for the island escarpment phenomenon. 
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CALCULATIONS [November 1978 Data] 


See the NOMENCLATURE section for a description of the symbols and units. 


AT1 = NQ)sSN(D6. NON gNO)9N(W) 19 


and 


N(1), Mode: 1 = Primary Sensor, 2 = Backup Sensor 
N(1)2 Resolution: 3 = N(1)7. N(1)g, 4= NO), . NO), 
N()); Interface Number: 1 = 1, 2=2, 3=3,4=4 
N(1)q Polarity; 1 = +, 2=— 

“N(1) 9 is rounded off to N(1), when N(1), =4 


AT2 = N(2)5N(2)¢ . N(2)7N(2)g.N(2)9N(2) 15 


and 


N(2), Mode: 1 = Primary Sensor, 2 = Backup Sensor. 
N(2)2 Resolution: 3 = N(2);.N(2)3, 4 = N(2),. N(2)7 
N(2); Interface Number: 1 = 1, 2=2, 3 =3, 4=4 
N(2)4 Polarity: 1 =+, 2=— 

*N (2) 19 is rounded off to N(2), when N(2), = 4 


If V(O5) < 2.570 volt, then: 
WS1 =—6.32834 x 1072 + 2.42269 x V(05) —0.399002 x V(05)2+8.72726 x 107? x V(05)? 


If V(05) > 2.570 volt, then: 
WS1 = 3.08914 x 1072+ 1.93793 x V(05) —6.72817 x 1074 x V(05)2 


If V(06) < 2.558 volt, then: 
WS2 = -7.13544 x 10-2 + 2.52039 x V(06) — 0.404911 x V(06)? + 8.22668 x 10-2 x V(06)? 


If V(06) > 2.558 volt, then: 
WS2 = 0.181091 + 1.92438 x V(06) + 1.55926 x 1074 x V(06)? 


v'(03) = V(03) + DPIFCAL 
DP1 =—20.825 + 6.25047 x V'(03) — 4.04968 x 107? x V'(03)? + 1.43719 x 107? x V'(03)° 
v'(04) = V(04) + DP2FCAL 
DP2 =—20.5158 + 6.06103 x V'(04) — 2.15130 x 10-2 x V'(04)? + 1.890459 x 107? x V'(04)? 
If |V(O1)| > 1.0 volt, then: 

bnew 
TUR) = |1.07722 x 10 


if |V(O1)| < 1.0 volt, then: 
TUR1 ="NO DATA" 


If | V(02)| > 1.0 volt, then: 








[- v(02 Il 
TUR2 = {1.07722 x 10' 


If |V(02)| < 1.0 volt, then: 
TUR2 = "NO DATA" 


If V(09) < 6.324 volt, then: 
WD = 145 + 34 x V(09) 


If V(O9) > 6.324 volt, then: 
WD = -215 + 34 x V(09) 


BP = 941.51 + 14.123 V(07) + 0.22644 x V(07)? — 1.41907 x 10-2 x v(07)3 
SR = 139.46 x V(08) 

WTB = 10 + V(10) 

Z1 = 18.450 

Z2 = 6.258 

ZB = 10.000 

ZTUR1 = Z1 + 0.60 

ZTUR2 = Z2 + 0.60 

GMH = VJZ1i x Z2 

ATK1 = AT1 + 273.160 

ATK2 = AT2 + 273.160 

DPK\ = DP1 + 273.160 

DPK2 = DP2 + 273.160 

Al = 373.160/DPK | 

A2 = 373.160/DPK2 
B\ = Ali — 1.000 
B2 = A2 — 1.000 | 
1 


A\ 





C1 = 1.000 — 








1 
ca 1000 - [5] 


GGC10 = —3.49149 
GGC11 = 11.344 
GGC12 = 5.02808 
GGC13 = —7.90298 
GGC14 = 8.1328 x 10 
GGC15 = -1.3816 x 107’ 
D1 = [10'66C1 (81) — 1 990 
¥ D2 = [10'66C1)*(82) — 1 999 
52 








| Fl= [19(GeCtD xc) — 1.000 
| F2 = [10°66C') (C2) _ 1 900 


VP =~ 1013.246 x (41) 662 x 1O!GGC 13x BI +IGGC 14x DIN +{GGC15x Fl} 
VP2 = 1013.246 x (A 2) GaC12 x 19! GGC 13x B2)+1GGC 14x D2)+1GGC15x F2] 
G1 = 373.160/ATK1 

G2 = 373.160/ATK2 

Q1 = G1 — 1.000 

Q2 = G2 — 1.000 


: 













1 
J1 = 1.000 [Z| 


a. 

G2 

Ki= [19°GEC10x(QI)} — 1.000 

K2 = [10°64'1*(27)) — 1.000 

M1 = (10°66C'DxUD) _ 1.900 

M2 = [10°66C') x22) _ 1.900 

SVP1 = 1013.246 x (G1) G62 x 10! GGC13x QI+IGGC14x K 1]+1GGC1Sx M1] 


J2 = 1.000 — | 


SVP2 = 1013.246 x (G2) 66C12 X JQ!GGC13x Q21+1GGC14* K 2} +1GGC1Sx M2} 
WVGC = 4.6150 x 10? 


























100 
re -pilicie Lwrattarey| 
100 
ee Lr teare | 
VPI 
RHI = |2F | x 100 
vP2 
RH? = | <7 x 100 
VPI 
sm =| Z| x 0.62197 
vP2 
SH? = |*5 | x 0.62197 





VTK1 = (ATK1) x [1.000 + (SH1 x 0.608)] 
VTK2 = (ATK2) x [1.000 + (SH2 x 0.608)]} 
VT1 = VTK1 — 273.160 
V72 = VTK2 — 273.160 


PTK1 = (ATK!) + (0.0098 x Z1) 
PTK2 = (ATK2) + (0.0098 x Z2) 
PT1 = PTK1 — 273.160 

PT2 = PTK2 — 273.160 

VPTK1 = (VTK1) + (0.0098 x Z1) 
VPTK2 = (VTK2) + (0.0098 x Z2) 
VPT1 = VPTK1 — 273.160 

VPT2 = VPTK2 — 273.160 


(WS1 — WS2) x Ln [23] 
* Po ge ee ener rr ees eee *See ERRATA NOTE 3, page 60. 
DWSDZ ul 
(SH1 — SH2) x Ln [24 
+ DSHDZ = ———__°\—* 
Zl 
(PTK1 — PTK2) x Ln [23| 
* DPTDZ = ————, 
Z\ 
(VPTK1 — VPTK2) x Ln | 20 
* DVPTDZ = oan 
pig co AE : ATK2 
ety se me : VPTK2 
GA = 9.7959 


(GA) x (DPTDZ) 


* RI = 
(VPTKM) x [DWSDZ]? 


If R7 < 0, then compute ZOL as a function of R/ solving the below equation in reverse via the 


Newton-Raphson Method for: 


{1 — 15(ZOL)]'? 


RI = (0.74(ZOL)] x 
e {1 — 9(ZOL)]'” 


Let J = R/as an arbitrary starting point and 


F(ZOL*) 


ZOL!*' = ZOL/ — 
F'(ZOL 


Where: 


ZOL! is the current guess. 
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ZOL/!*' is the next iteration. 


and 


74(ZOL){1 ~ 1S(ZOL)\'?___p, 
[1 — 9(ZOL)]'? 


F(ZOL) = [i741 15(Z0L)}")| - eee | + faceonae | 


F(ZOL) = 


[2-30(ZOL)]"”” (2-18(ZOL)]'? 
Execute the iterations until: 


J+l J 
ZOL ZOL < 0005 


|ZOL/*!! 


At which point an exceptable value for ZOL has been computed. 





Then: 
GMH 
MOL ZOL 
a ZI 
Z\OL = MOL 
Z2 
Z20L MOL 


W1 = [1 — 15(Z10L)]"4 
W2 = [1 — 15(Z20L)]"4 
Y1 = [1 —9(Z10L)]'” 
¥2 = [1 — 9(Z20L)]'”” 


With arctan expressed in radians: 


























PSIVZ1 = 2x Ln) *+"1) + tn Latte — (2 arctan(W1)) + = 
PSI1Z2 =2 x Ln im +Ln Ler — [2 arctan(W2)] + = 
PSI2Z1 = Ln —_— 
PSI2Z2 = Ln in 








If RI > 0, then compute ZOL as a function of R/ solving the below equation in reverse via the 


Quadratic Solution for: 





_ 9.74(ZOL) + 4.1(ZOL)? 


id {1 + 4.7(ZOL)}* 


Let: 


—B —~VB?-—4(A x C) 
201 2A) 


Where: 


A = (22.09 x RI) — 4.7 
B = (9.4 x RI) — 0.74 
C=RI 

Then: 


GMH 
ZOL 
Zi 
MOL 
ae. 
MOL 
PSI1Z1 = —4.7(Z10OL) 


PSI1Z2 = —4.7(Z20L) 
—4.7(Z10L) 
0.74 
—4.7(Z20L) 
0.74 


MOL = 
Z10OL = 


Z20L = 


PSI2Z1 = 


PSI2Z2 = 


ZB 
ZOLT MOL 


VKC = 0.35 


If (WS1 — WS2) =0, then set (WS1 — WS2) = 0.001 


([Ln(Z1) — PSI1Z1) — [Ln(Z2) — PSI1Z2) 


ae (WS1 — WS2) 


If (SH1 — SH2) = 0,then set (SH1 — SH2) =1 x 107’ 


[Ln(Z1) — PSI2Z1) — [Ln(Z2) — PSI2Z2) 
(SH1 — SH2) 


If (PTK1 — PTK2) =0, then set (PTK1 — PTK2) = 0.0001 


{[Ln(Z1) — PS12Z1) — [Ln(Z2) — PSI2Z2) 
(PTK1 — PTK2) 


If (VPTK1 — VPTK2) = 0, then set (VPTK1 — VPTK2) = 0.0001 


[Ln(Z1) — PS12Z1) — [Ln(Z2) — PSI2Z2) 
(VPTK1 — VPTK2) 


If [Ln(TUR1) — Ln(TUR2)] = 0, then set [Ln(TUR1) -- Ln(TUR2)} 
= Ln(1 x 1075) = -11.51292546 


SHSLOP = 


PTKSLOP = 


VPTKSLOP = 








Pee ee ee 


Te a ee a 





[Ln(ZTUR1) — Ln(ZTUR2)]) 
[Ln(TUR1) — Ln(TUR2)} 
VKC 
WSSLOP 


Fricv) 
ZB 


TURSLOP = 


FRICV = 


DRAG = | 





WS1 + WS2 


W: = 
‘SM 


_ | (WSM) x (VKC) 
RL = Ln(GMH) | ae | 


If RL < —227.95, then: 
RUFL =} x 10-% 
If RL > —227.95, then: 


RUFL = exp(RL) 
VKC 
0.74(SHSLOP) 
VKC 
0.74(PTKSLOP) 
VKC 
0.74( VPTKSLOP) 
VTK1 + VTK2 
2 
.34838(BP) 
VTKM 


SCLSH = 
SCLPT = 
SCLVPT = 
VTKM = 


AIRD = 


SPHEAT = 0.240|1 + 0.90 x 10° 





a 


MOFLX = —(AIRD) x (FRICV)? 

HUFLX = —(AIRD) x(FRICV) x (SCLSH) 

HEFLX = —(AIRD) x (FRICV) x(SCLPT) x {(SPHEAT) x 4.18684) 
ATM = ATKM — 273.160 

LHV = {597.31 — [.56525(ATM)]} x 10° 

BOYFLX = HELFX + ((HUFLX) x (LHV) x (4.18684)] 

THFLX = BOYFLX — SR 

LAMBDA = 5.4 x 10~® 


ie 
OPTIR1 = - (nse «fs 7.53 x10°_ SEIS | x crurn x 10- 


(ATK 1)? (LAMBDA)? 
71.6(BP) 1.83 «10° |’ aa 
OPTIR2 = TUR2) x 10 
si [scars SOP) «i+ (LAMBDA) | eer 


WSB = WS)| — [.61195 (WS1 — WS2)) 
DPKB = DPK\ — [.53098(DPK1 — DPK2)] 
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ATKB = ATK) — [.55743 (ATK1 — ATK2)] 


TUR | 
TUR2 


TRT = Ln (TUR}1) - sso x Ln | 


TURT = exp(TRT) 

DPB = DPKB — 273.160 
ATB = ATKB — 273.160 

AB = 373.160/DPKB 

BB = AB — 1.000 

ak 

AB 

DB = [10(6C10(88)) _ 1.990 


FB = (10(GGC11)x(CB)) — 1.000 
VPB = 1013.246 x (AB) G6612 x 1Q!GGC13x BB)+1GGC 14x DB) +(GGC15x FB] 


CB = 1.000 — | 


GB = 373.160/ATKB 
QB = GB — 1.000 


1 
JB = 1.000 [35 


KB = [10'GEC 10) x(08)) — 1,000 
MB = [10°GEC1D x8) — 1.000 
SVPB = 1013.246 x (GB) &60!2 x 10! GGC13x QB)+1GGC14x KB] +[GGC15x MB] 


100 


AHB = (VPB) X | GC) x (ATKB) 


VPB 
SVPB 


VPB 
SHB Bp | * 0.62197 


VTKB = (ATKB) x [1.000 + (SHB x 0.608)] 
VTB = VTKB — 273.160 

PTKB = (ATKB) + [0.0098 x ZB] 
PTB = PTKB — 273.160 

VPTKB = (VTKB) + [0.0098 x ZB] 
VPTB = VPTKB — 273.160 

WTKB = WTB + 273.160 

ATWTB = ATKB — WTKB 

PTWTB = PTKB — WTKB 

VTWTB = VTKB — WTKB 
VPTWTB = VPTKB — WTKB 


RHB = «100 











SHTC = 0.92 x 1073 


SPHEATB = 0.240[1 + 0.90 (SHB)} x 10° 


34838 (BP) 
VTKB 


VVPTC = 0.002 — [(SHTC) x (WSB) x (PTWTB)] 
HEFLXB = (AIRDB) x (VVPTC) x ((SPHEATB) x 4.18684] 
AW = 373.160/WTKB 

BW = AW — 1.000 


AIRDB = 


1 
cw = 1.000 — [1 


DW = [10°G6C10)x(aW} — 1.000 


FW = [10(Gectnx(cw) — 1.000 
VPW = 1013.246 x (A W)GGC12 x 1Q'GGC13x BW) +1GGC14x DW) +1GGC15x Fw] 


AHW = (VPW) + aac arrear| 

MTC = 1.32 x 1073 

VVHC = (MTC) x (WSB) x [AHW — AHB] 

HUFLXB = VVHC 

LHVB = (597.31 — [.56525(ATB)]} x 10° 

BOYFLXB = HEFLXB + ((HUFLXB) x (LHVB) x (4.18684)} 
THFLXB = BOYFLXB — SR 

DRAGB = {0.63 + [0.066( WSB)]} x 10-3 

RSTRES = (DRAGB) x (WSB)? x (AIRDB) 

VVLC = RSTRES/AIRDB 


MOFLXB = —(RSTRES) 


rucre = « EEE 

RLB = Ln (zB) ~ <¥SB) AKO) 
If RLB < -227.95, then: 

RUFLB = 1 x 10- 

If RLB > -227.95, then: 

RUFLB = exp(RLB) 


—HUFLXB 
(AIRDB) x (FRICVB) 
—HEFLXB 
(AIRDB) x (FRICVB) x [4.18684(SPHEATB)} 


SCLSHB = 


SCLPTB = 
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, 
' 








~BOYFLXB \ 


*“SCLVPTB = "See ERRATA NOTE 2 on page 60. 


(AIRDB) x(FRICVB) x [4.18684(SPHEATB)] 


—(ATKB) x [FRICVB)} 


MOLB = oA) x (VVPTC) x(VKO) 
GMH 
ZOLB = tO 


If ZOLB < 0, then: 


{1 — 15(ZOLB)]'” 


RIB = '(0.TZOLEN gar a 


If ZOLB > 0, then: 
Si) (ZOLB) x [0.74 + 4.7(ZOLB)} 


RIB 
[1 + 4.7(ZOLB)}? 

ZB 

20LTR < = 
2 2 

(77.6) x (BP) 7.53 x 10° = 

OPTIRT = |-—-————|_ x |1 + —=———— | x (TuRT) x 10 
(ATKB)? | 7s | ; ; 








ERRATA NOTES 


1 The Bulk Water Temperature reading of 15.667°C for Run Number 7811140800 is incorrect. It should have 
read 15.094°C. 


z The equation for "SCLVPTB" above is incorrect, since the definition for "BOYFLXB" had been redefined. The 
value presented for the bulk derived Scaling Virtual Potential Temperature is incorrect. 


3. The general form of the approximation for the partial derivatives employed in computing the profile derived 
Richardson number stability were obtained from a key article by C. A. Paulson, “The Mathematical Representa- 
tion of Wind Speed and Temperature Profiles in the Unstable Atmospheric Surface Layer,” J. of Appl. Meteor., 
Vol. 9, pp. 857-861. Dec. 1970. The Paulson paper contains a typographical error in the equation for the partial 
derivative which was subsequently transcribed with the error into our calculations. The term Ln (Z1/Z2), in 
our subscript notation, which appears in the numerator of the partial derivative equations, should have appeared 
in the denominator. The value of Li (Z1/Z2) is a constant equal to 1.0812037. The correct values for the vari- 
ous partial derivatives (DWSDZ. DSHDZ, DPTDZ. and DVPTDZ) may be obtained by multiplying the 
presented value by 0.85543. The correct value for the profile Richardson number (RI) may be obtained by mul- 
lplying the presented value by 1.1690. These corrections have little practical effect on the conclusions arrived at 
in this report 
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NOMENCLATURE 
A = 

Al - 

A2 - 

AB - 
AIRD kg/m? 
AIRDB kg/m? 
AHI kg/m? 
AH2 kg/m? 
AHB kg/m? 
AHW kg/m? 
arctan( ) radians 
ATI . 
AT2 ‘c 
ATB °C 
ATK1 °K 
ATK2 °K 
ATKB °K 
ATKM °K 
ATM *€ 
ATWTB °K 
AW - 

B am 





{November 1978 Data] 


Description: © 


Variable in quadratic solution for profile derived Z/L (or f) stability 
parameter when Richardson’s number > 0. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z1 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z2 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the 10 meter height. 


Mean air density at the Z1 and Z2 heights. 

Air density at 10 meter height. 

Absolute humidity at the Z1 height. 

Absolute humidity at the Z2 height. 

Absolute humidity at the 10 meter height. 

Absolute humidity at water level. 

Arctangent of (). 

Air temperature in °C at the Z1 height. 

Air temperature in °C at the Z2 height. 

Air temperature in °C at the 10 meter height. 

Air temperature in °K at the Z1 height. 

Air temperature in °K at the Z2 height. 

Air temperature in °K at the 10 meter height. 

Mean air temperature in °K of the measurements at the Z1 and Z2 heights. 
Mean air temperature in °C of the measurements at the Z1 and Z2 heights. 
Air temperature at the 10 meter height minus the bulk water temperature. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at water level. 


Variable in quadratic solution for profile derived Z/L (or %) stability 
parameter when Richardson’s number > 0. 
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Bl 


B2 


BB 


BOYFLX 


BOYFLXB 


BP 


BW 


Cl 


C2 


CB 


CW 


D1 


D2 


DB 


DPI 
DPIFCAL 
DP2 
DP2FCAL 


DPB 


watt/m? 


watt/m? 


mb 


‘¢ 
volt 
a % 
volt 


°C 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z1 meter height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z2 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the 10 meter height. 


Profile derived sensible and latent heat flux. This is nor 
buoyancy flux, for which the symbol was previously employed. 


Bulk aerodynamic derived sensible and latent heat flux. This is nor 
buoyancy flux, for which the symbol was previously employed. 


Mean barometric pressure at the Z1, Z2, and 10 meter heights. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at water level. 


Variable in quadratic solution for profile derived Z/L (or &) stability 
parameter when Richardson’s number > 0. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z1 height. 





Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z2 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the 10 meter height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the water level. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z1 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z2 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the 10 meter height. 


Dew point in °C at the Z1 height. 
Dew point field calibration correction at the Z1 height. 
Dew point in °C at the Z2 height. 
Dew point field calibration correction at the Z2 height. 


Dew point in °C at 10 meter height. 
62 


en ee ee -—- 





DPK1 
DPK2 
DPKB 


DPTDZ 


DRAG 


DRAGB 


DSHDZ 


DVPTDZ 


DW 


DWSDZ 


exp() 


FI 


F2 


FB 


FRICV 
FRICVB 


FW 


Gl 





°K 
°K 
°K 


°K/m 


kg/kg -m 


°K/m 


m/sec - m 


m/sec 


m/sec 


Dew point in °K at the ZI height. 
Dew point in °K at the Z2 height. 
Dew point in °K at the 10 meter height. 


Profile derived partial derivative of potential temperature 
with respect to height. 


Profile derived atmospheric drag coefficient at the 
10 meter height. 


Bulk aerodynamic derived atmospheric drag coefficient 
at the 10 meter height. 


Profile derived partial derivative of specific humidity 
with respect to height. 


Profile derived partial derivative of virtual potential 
temperature with respect to height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at water level. 


Profile derived partial derivative of wind speed with respect 
to height. 


Exponent to base of natural logarithm e of ( ), e = 2.71828... 





Variable in Newton-Raphson method solution for profile derived Z/L (or £) 
stability parameter when Richardson’s number < 0. 


Variable in Newton-Raphson method solution for profile derived Z/L (or 2) 
stability parameter when Richardson’s number < 0. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z1 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the Z2 height. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at the 10 meter height. 


Profile derived friction (or scaling) velocity. 
Bulk aerodynamic derived friction (or scaling) velocity. 


Intermediate step in Goff-Gratch formulation of vapor pressure 
at water level. 


Intermediate step in Goff-Gratch formulation 
of saturated vapor pressure at the Z1 height. 
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G2 


GA 


GB 


GGC10 
GGCll 
GGCl12 
GGC13 
GGC14 
GGC15 


GMH 


HEFLX 
HEFLXB 
HUFLX 
HUFLXB 


J 


Ji 


J2 


JB 


Kl 


K2 


KB 


m/sec? 


meter 


watt/m2 
watt/m? 
kg/sec - m? 


kg/sec - m? 





Intermediate step in Goff-Gratch formulation 
of saturated vapor pressure at the Z2 height. 


Gravitational acceleration constant; set equal to 
9.7959 for the San Nicolas Island latitude. 


Intermediate step in Goff-Gratch formulation 
of saturated vapor pressure at the 10 meter height. 


Goff-Gratch formulation constant: —3.49149. 
Goff-Gratch formulation constant: 11.344. 
Goff-Gratch formulation constant: 5.02808. 
Goff-Gratch formulation constant: —7.90298. 
Goff-Gratch formulation constant: 8.1328 x 107°. 
Goff-Gratch formulation constant: —1.3816 x 1077. 
Geometric mean of heights Z1 and Z2. This is the 
equivalent to the mean measurement height as 
plotted on a Ln (height) vs. linear graph. 

Profile derived sensible heat flux. 

Bulk aerodynamic derived sensible heat flux. 
Profile derived specific humidity flux. 


Bulk aerodynamic derived specific humidity flux. 


Variable in Newton-Raphson method solution for profile derived Z/L (or ¢) 
Stability parameter when Richardson's number < 0. 


Interm ‘iate step in Goff-Gratch formulation of saturated vapor 
pressure at the Z1 height. 


Intermediate step in Goff-Gratch formulation of saturated vapor 
pressure at the Z2 height. 


Intermediate step in Goff-Gratch formulation of saturated vapor 
at the 10 meter height. 


Intermediate step in Goff-Gratch formulation of saturated vapor 
pressure at the Z1 height. 


Intermediate step in Goff-Gratch formulation of saturated vapor 
pressure at the Z2 height. 


Intermediate step in Goff-Gratch formulation of saturated vapor 
pressure at the 10 meter height. 
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FO Se eee ee 


LAMBDA 


LHV 
LHVB 
Ln() 


M1 


M2 


MB 


MOFLX 
MOFLXB 
MOL 
MOLB 


MTC 


N(), 


OPTIR1 


OPTIR2 


OPTIRT 


PSI1Z1 


PSI1Z2 


PSI2Z1 


PS1I2Z2 


PTI 


meter 


ITcal/kg 


ITcal/kg 


Nt/m? 
Nt/m?2 


meter 


°K/m-23 


°K/m-23 


°K/m-23 


Wavelength employed for computing optical refractive index; 
set equal to 5.4 x 107°. 


Profile derived latent heat of water vaporation. 
Bulk aerodynamic derived latent heat of water vaporation. 
Natural logarithm of ( ) with base e = 2.71828 .... 


Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the Z1 height. 


Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the Z2 height. 


Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the 10 meter height. 


Profile derived momentum flux. 

Bulk aerodynamic derived momentum flux. 
Profile derived Monin-Obukhov length. 

Bulk aerodynamic derived Monin-Obukhov length. 


Bulk aerodynamic moisture transfer coefficient; 
set equal to 1.32 x 107°. 


Digital channel number () consisting of 10 parallel digital characters, 
n=1,2,...10. 


Optical refractive index parameter (C7) 
at the ZTURI height. 


Optical refractive index parameter (C,?) 
at the ZTUR2 height. 


Optical refractive index parameter (C,2) 
at the 10 meter height. 


Businger wind speed profile stability correction 


(W,) for the Ln (Z1) height. 


Businger wind speed profile stability correction 
(W,) for the Ln (Z2) height. 


Businger temperature and humidity stability correction 
(W,) for the Ln (Z1) height. 


Businger temperature and humidity stability correction 
(W,) for the Ln (Z2) height. 


Potential temperature in °C at the Z] height. 
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PT2 °€ Potential temperature in °C at the Z2 height. 

PTB °C Potential temperature in °C at the 10 meter height. 

PTK1 °K Potential temperature in °K at the Z1 height. 

PTK2 °K Potential temperature in °K at the Z2 height. 

PTKB °K Potential temperature in °K at the 10 meter height. 
PTKSLOP - Profile derived physical slope of Ln (height) with Businger 


stability corrections vs. potential temperature. 


PTWTB °K Potential temperature at the 10 meter height 
minus the bulk water temperature. 


Ql - Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the Z1 height. 


Q2 - Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the Z2 height. 


QB - Intermediate step in Goff-Gratch formulation of saturated 
vapor pressure at the 10 meter height. 


RH1 % Relative humidity at the Z1 height. 

RH2 % Relative humidity at the Z2 height. 

RHB % Relative humidity at the 10 meter height. — 

RI - Profile derived gradient Richardson number stability. 

RIB - Bulk aerodynamic derived gradient Richardson number stability. 
RL - Intermediate step in computation of profile derived 


roughness length. 


RLB - Intermediate step in computation of bulk aerodynamic 
derived roughness length. 


RSTRESS kg/m - sec? Bulk aerodynamic derived Reynold’s stress. 


RUFL meter Profile derived roughness length. 

RUFLB meter Bulk aerodynamic derived roughness length. 

SCLPT °K Profile derived scaling potential temperature. 

SCLPTB °K Bulk aerodynamic derived scaling potential temperature. 
SCLSH kg/kg Profile derived scaling specific humidity. 

SCLSHB kg/kg Bulk aerodynamic derived scaling specific humidity. 
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SCLVPT 
SCLVPTB 
SHI 

SH2 

SHB 


SHSLOP 


SHTC 


SPHEAT 
SPHEATB 
SR 

SVP1 
SVP2 
SVPB 
THFLX 
THFLXB 


TRT 


TTABLE 


TURI 


TUR2 


TURSLOP 


TURT 


v() 


vi () 


°K 
°K 
kg/kg 
kg/kg 


kg/kg 


ITcal/kg °K 
ITcal/kg °K 
watt/m? 

mb 

mb 

mb 
watt/m? 


watt/m2 


meter 


°K/m-23 


°K/m-23 


°K /m-2/3 


volt 


volt 


Profile derived scaling virtual potential temperature. 

Bulk aerodynamic derived scaling virtual potential temperature. 
Specific hunnidiey at the Z1 height. 

Specific humidity at the Z2 height. 

Specific humidity at the 10 meter height. 


Profile derived physical slope of Ln (height) with Businger 
stability corrections vs. specific humidity. 


Bulk aerodynamic sensible heat transfer coefficient; 
set equal to 0.92 x 107°. 


Profile derived specific heat of air. 

Bulk aerodynamic derived specific heat of air. 
Sky and solar radiation flux. 

Saturated vapor pressure at the Z1 height. 
Saturated vapor pressure at the Z2 height. 
Saturated vapor pressure at the 10 meter height. 
Profile derived total heat budget flux. 

Bulk aerodynamic derived total heat budget flux. 
Intermediate step in computing the microthermal 
turbulence temperature structure parameter (C,’) 
at the 10 meter height from profile C,? measurements. 


Tide table referenced to mean sea level. 


Microthermal turbulence temperature structure parameter 
(C?) at the ZTURI height. 


Microthermal turbulence temperature structure parameter 
(C?) at the ZTUR2 height. 


Profile derived slope of Ln (height) vs. Ln (microthermal turbulence 
temperature structure parameter (C,’)). 


Microthermal turbulence temperature structure parameter 
(C?) at the 10 meters height. 


Average raw data voltage for analog channel number ( ). 


Field calibration corrected average data voltage for analog 
channel number ( ). 
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VKC - Von Kaérmdn constant; set equal to 0.35. 

VPI mb Vapor pressure at the Z1 height. 

VP2 mb Vapor pressure at the Z2 height. 

VPB mb Vapor pressure at the 10 meter height. 

VPT1 %€ Virtual potential temperature in °C at the Z1! height. 

VPT2 fb Virtual potential temperature in °C at the Z2 height. 
VPTB °C Virtual potential temperature in °C at the 10 meter height. 

VPTK1 °K Virtual potential temperature in °K at the Z] height. 

VPTK2 °K Virtual potential temperature in °K at the Z2 height. 

VPTKB °K Virtual potential temperature in °K at the 10 meter height. 

VPTKM °K Mean virtual potential temperature of the measurements 


at the Z1 and Z2 heights. 


VPTKSLOP — Profile derived phsycial slope of Ln (height) with Businger 
Stability corrections vs. virtual potential temperature. 

VPTWTB °K Virtual potential temperature at the 10 meters height 
minus the bulk water temperature. 

VPW mb Vapor pressure at water level. 

VTl °C Virtual temperature in °C at the Z1 height. 

VT2 be @: Virtual temperature in °C at the Z2 height. 

VTB °© Virtual temperature in °C at the 10 meter height. 

VTK1 °K Virtual temperature in °K at the Z1 height. 

VTK2 °K Virtual temperature in °K at the Z2 height. 

VTKB °K Virtual temperature in °K at the 10 meter height. 

VTKM °K Mean virtual temperature of measurements at Z1 and Z2 heights. 

VTWTB °K Virtual temperature at the 10 meter height minus 


the bulk water temperature. 


VVHC m - kg/sec: m? Bulk aerodynamic derived vertical velocity and absolute 
humidity covariance. 


VVLC m?/sec? Bulk aerodynamic derived vertical velocity and longitudinal 
velocity covariance. 
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VVPTC 


WD 
ws! 
Wws2 
WSB 
WSM 


WSSLOP 


WTB 
WTKB 
WvVGC 


wi 


w2 


Yl 


Y2 


ZI 


Z2 


ZB 


Z1OL 
Z20L 
ZOL 


m °K/sec 


degree 
m/sec 
m/sec 
m/sec 


m/sec 


*€ 
°K 


Nt - m/kg °K 


meter 


meter 


meter 


Bulk aerodynamic derived vertical velocity and potential 
temperature covariance. 


Wind direction referenced to true North. 

Wind speed at the Z1 height. 

Wind speed at the Z2 height. 

Wind speed at the 10 meter height. 

Mean wind speed of the measurements at the Z1 and Z2 heights. 


Profile derived physical slope of Ln (height) with Businger 
stability corrections vs. wind speed. 


Bulk water temperature in °C. 
Bulk water temperature in °K. 
Water vapor gas constant; set equal to 4.6150 x 107. 


Intermediate step in computation of the Businger 
stability correction for wind speed profile when 
Richardson number < 0 at the Z1 height. 


Intermediate step in computation of the Businger 
Stability correction for wind speed profile when 
Richardson number < 0 at the Z2 height. 


Intermediate step in computation of Businger stability 
correction for temperature and humidity profiles 
when Richardson number < 0 at the Z1 height. 


Intermediate step in computation of Businger stability 
correction for temperature and humidity profiles 
when Richardson number < 0 at the Z2 height. 


Upper level average measurement height of profile instruments 
(except C,?) referenced to ground surface. 


Lower level average measurement height of profile instruments 
(except C,?) referenced to ground surface. 


This is the standard height for the bulk aerodynamic calculations; 
set equal to 10 meters and referenced to ground surface. 


Profile derived Z/L (or {) stability parameter computed for the Z1 height. 


Profile derived Z/L (or &) stability parameter computed for the Z2 height. 


Profile derived Z/L (or £) stability parameter computed 
for the geometric mean height which is the relevant height for the 
conversion of Richardson number to Z/L. 
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ZOLB 


ZOLT 


ZOLTB 


ZTURI 


ZTUR2 


P Any 


meter 


meter 


Bulk aerodynamic derived Z/L (or {) stability 
parameters computed for the geometric mean height. 


Profile derived Z/L (or {) stability 
parameter computed for the 10 meter height. 


Bulk aerodynamic derived Z/L (or {) stability 
parameter computed for the 10 meter height. 


Upper level measurement height of microthermal 
turbulence temperature structure (C,2) sensors 
referenced to the ground surface. 


Lower level measurement height of microthermal 
turbulence temperature structure (C,2) sensors 
referenced to the ground surface. 
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COMPUTER PROGRAM VARIABLE AND DATA ARRAY LABELS 


See the NOMENCLATURE section for a definition of the parameter symbols. 


I. "CRNCH2" PROGRAM: 
Data Array Label: 


D(1) 
D(2) 
D(3) 
D(4) 
D(5) 
D(6) 
D(7) 
D(9) 
D(10) 
D(i2) 
D(18) 
D(19) 


Variable Label: 


ri 

r2 

3 

r4 

r5 

r6 

7 

8 

r9 

rl0 
rll 
ri2 
r13 
rl4 
r16 
rl7 
r18 
ri9 
r20 
r21 
r22 
123 
r24 
125 
126 
r27 
r28 
r29 
r30 
31 





Parameter: 


TURI 
TUR2 
DP1 
DP2 
WwS1 
WSs2 
BP ' 
WD 
WTB 
SR 
AT} 
AT2 


Parameter: 


ATI 
AT2 
OP! 


VP1 
VP2 
SVP1 
SVP2 
AHI 
AH2 
RH1 
RH2 
SH] 
SH2 
OPTIRI 
OPTIR2 
VTl 
VT2 
PTi 
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| 
| 
| 
| 
| 


132 PT2 











133 VPT1 

134 VPT2 

135 DWSDZ 
136 DSHDZ 
137 DPTDZ 
138 DVPTDZ 
r40 RI 

r41 ZOL 

r42 MOL 

143 Z1OL 

144 Z20L 

r45 PSIIZ1 

146 PSI1Z2 
147 PSI2Z1 

r48 PS12Z2 
149 FRICV 
r50 RUFL 

r51 SCLSH 
r52 SCLPT 
153 SCLVPT 
154 AIRD 

r55 MOFLX 
156 HUFLX 
r57 HEFLX 
158 THFLX 
r59 WSSLOP 
r60 SHSLOP | 
r61 PTKSLOP | 
162 VPTKSLOP | 
163 TURSLOP 
164 GMH 

165 Z) 

166 Z2 

169 DRAG 
r70 ZOLT 

71 SPHEAT 





II. "CRNCH3" PROGRAM 





Data Array Label: Parameter: 
D(10) WTB 
F | Variable Label: Parameter: 
9 WTB 
r12 SR 
rl3 ATB 
q rl4 WSB 
ri5 DPB 
| rl6 TURT 
\ rl8 ATWB 


| r19 ' PTWTB 
| 490 








| 
] 
| 
/ 


120 VTWTB 
121 VPTWTB 
r22 ZB | 
123 PTB | 
124 VTB ] 
r25 VPTB 
126 AHB | 
127 RHB 
128 SHB 
129 VPB | 
130 SVPB 
131 OPTIRT 
135 RSTRES i 
136 VVHC | 
138 VVPTC 
r40 RIB ] 
r41 ZOLB | 
r42 MOLB | 
r49 FRICVB 
r50 RUFLB 
51 SCLSHB 
r52 SCLPTB 
r53 SCLVPTB | 
154 AIRDB | 
155 MOFLXB 
156 HUFLXB | 
157 HEFLXB 
158 THFLXB | 
160 LVHB | 
r61 VPW 
162 ; AHW 
: 164 GMH 
r69 DRAGB 
r70 ZOLTB { 
r71 SPHEATB 
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